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THEORETICAL  PREDICTK^ 

OF  VIBRATIONAL  CIRCULAR  DIOIROISM  SPECTRA 


1.  INTRODUCTION 

A  requirement  for  this  study  was  the  interest  of  the  Detection  Group  at  CRDEC  in  being  able  to 
detect  biological  materials  in  the  field.  A  reasonable  means  to  this  objective  is  through  the 
recognition  that  sugars  arc  distinguishing  features  of  biological  materials.  Sugars  arc  optically  active 
and  can  be  studied  through  a  variety  of  techniques.  The  technique  of  interest  here  is  vibrational 
circular  dichroism  (VCD).  The  predicted  VCD  spectra  can  also  be  related  to  the  (1,4)  matrix  element 
of  the  Mueller  matrix,  to  be  discussed  in  section  II.  A  procedure  to  obtain  all  of  16  elements  of  the 
(4x4)  Mueller  matrix  simultaneously  by  a  two  laser  spectroscopic  measurement  is  presently  being 
implemented  by  one  of  the  authors  (AHC).  The  calculated  spectra  of  sugars,  at  a  minimum,  should 
assist  laboratory  measurements  on  sugars.  The  practical  questions  as  to  whether  real  biological 
materials  can  be  detected  and  the  limits  of  detectability  await  future  experimental  studies.  Because 
the  techniques  of  computational  chemistry  are  now  cheaper  and  faster  than  extensive  experimental 
techniques,  quantum  chemistry  is  employed  to  identify  spectral  regions  of  interest  and  what 
compounds  can  be  detected  by  VCD.  Furthermore,  a  library  can  be  constructed  that  will  permit 
matching  of  experimental  results  to  unambiguously  identify  unknown  compounds.  This  report  will 
focus  on  the  application  of  current  computational  techniques  to  the  study  of  sugars,  sugar-like 
molecules  and  other  molecules  of  interest.  This  study  can  be  viewed  as  the  beginning  phase  of  the 
generation  of  a  library  of  theoretically  predicted  VCD  spectra. 

This  report  is  organized  in  the  following  sections:  section  II  is  a  background  section  on  VCD  in 
general,  the  Stephens’  formulation  (1,2)  of  the  rotational  strength  of  a  vibrational  transition  and  the 
concept  of  a  Mueller  matrix;  section  III  focusses  on  three  classes  of  molecules  that  were  investigated; 
section  IV  discusses  the  work  suggested  for  continuation  of  this  study;  section  V  summarizes  the 
report. 


2.  BACKGROUND 


2.1  Vibrational  Circular  Dichroism 

Molecules  that  interact  differently  with  left  circularly  polarized  (LCP)  and  right  circularly 
polarized  (RCP)  radiation  arc  optically  active.  Two  phenomena  linked  to  optical  activity  arc: 

( 1 )  optical  rotation 

(2)  circular  dichroism 

Optical  rotation  involves  the  rotation  of  a  plane  of  linearly  polarized  radiation  as  it  is  passed  through  a 
sample.  Circular  dichroism  is  the  difference  in  the  absorption  of  LCP  and  RCP  radiation  by  a  sample. 
In  the  absence  of  applied  fields,  optical  activity  arises  from  chirality,  which  exists  when  a  molecule  is 
not  superimposable  on  its  mirror  image  -  i.e.  a  molecule  must  not  possess  any  elements  of  symmetry  to 
be  optically  active. 

In  the  early  1970’s  instrumentation  was  devised  to  measure  the  VCD  spectrum  of  a  chiral 
molecule.  Here,  the  difference  in  absorptivity  of  a  sample  between  LCP  and  RCP  radiation  is  plotted 
as  a  function  of  wavenumber.  A  typical  VCD  spectrum  is  shown  in  Figure  1.  Figure  1  includes  both 
the  VCD  spectrum  of  methyl  thiirane  over  the  wavenumber  range  of  1500-800  cm  *  and  the  IR 
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spectrum  for  the  purpose  of  comparison.  As  can  be  noted  the  rotational  strengths,  i.e.  the  difference 
in  absorptivity  of  a  sample  between  LCP  and  RCP  radiation,  can  be  positive  or  negative  and  the 
intensities  are  not  related  to  the  infrared  and  Raman  intensities.  The  essential  points  to  be  studied  in 
this  investigation  are  frequencies,  magnitudes  of  intensities  and  signs  of  intensities.  All  notations  used 
in  the  following  sections  are  consistent  with  those  of  Stephens  (1,2). 

2.2  Stephens’  Formulation  of  the  Rotational  Strength 

P.  J.  Stephens  (1,2)  has  developed  a  theory  for  the  determination  of  the  rotational  strength  of  a 
vibrational  transition  g-+  e,  R(g— ►  e),  for  a  chiral  molecule; 

R(^-^)=E  l/3=x,y,z]  (1) 

where  g  and  e  are  ground  and  excited  states  and  and  u  ^  are  electric  and  magnetic  dipole  moment 
operators.  For  a  given  fundai.iental  vibrational  mode  of  frequency  Wj,  the  electric  and  magnetic  dipole 
transition  moments  for  the  0-»  1  vibrational  transition  is  given  by; 


<o\(p^)l,\l>Mh  5xa.  (2) 

Aa 

<o\{p.„^)0\  1  >,=-(2h='w.)i/=E:  5,,^,  (3) 

Aa 

P  ^and  M  ^are  atomic  polar  and  atomic  axial  tensors  defined  by 

(4) 

(5) 

E^asHid/dXxa)  (6) 

(7) 

1^00  =  <dUR)/dX^J-^  I  {d^dM„H0)ldH0)u^o>  (8) 

l^hc)Y,^^Zye)R%  (9) 

7 


where 


R 


the  nuclear  geometry 

the  equilibrium  nuclear  geometry 

cartesian  displacement  coordinate'(a  =  x,y,z)  from  equilibrium 
of  nucleus  A 

charge  on  nucleus  A 
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B  equilibrium  position  of  nucleus  A 

G 

B  electronic  ground  state 

m:, 

B  electronic  contribution  to 

=  electronic  wavefunction  of  G  at  in  the 

presence  of  the  perturbation  H  b  -(/i„^*)^ 

B  antisymmetric  unit  third  rank  tensor 

Note  that; 


^  Xa—^S  Xa,tQi 
i 

i.e.  Sj^  is  a  matrix  element  in  the  transformation  between  the  cartesian  coordinates  and  the  normal 
coordinates. 


2.S  Mueller  Matrix 

Let  us  now  turn  to  a  brief  discussion  of  the  definition  and  use  of  a  Mueller  matrix.  First  let  us 
consider  the  transmission  of  electromagnetic  radiation.  The  interaction  of  polarized  radiation  with  an 
optical  element  (polarizer,  retarder,  reflector,  scatterer)  can  be  described  by  specifying  the 
polarization  state  of  incident  and  transmitted  radiation.  A  typical  experimental  setup  is  shown  in 
Figure  2.  One  convenient  description  of  the  radiation  involves  use  of  Stokes  vectors  and  the  use  of  a 
transformation  between  the  transmitted  radiation,  and  the  incident  radiation  a,  which  is  termed  the 
Mueller  matrix  (3,4),  T 


<T,=Tcr,- 


(10) 


All  of  the  information  concerning  the  transmission  of  radiation  by  an  optical  element  is  contained  in 
the  16  elements  of  the  transformation  Mueller  matrix.  The  Stokes  parameters,  components  of  a  Stokes 


vector  are  defined  by 

/=<£  £*  +£  £*  >  (11) 

Q=<E  £*  -£  £*  >  (12) 

U=<E  £'+££*>  (13) 

V=i<E  £*  -£  E*  >  (14) 


where  the  angular  brackets  indicates  a  time  average  over  an  interval  long  compared  with  the  period  of 
a  quasi-monochromatic  beam  described  by  an  electric  field  £  with  orthogonal  axes  perpendicular  to 
the  direction  of  propagation.  Typical  Stokes  vectors  arc  shown  in  Table  1. 

Now  if  scattered  radiation  from  a  single  particle  or  a  collection  of  particles  is  studied  as  shown  in 
Figure  3  an  analogous  Mueller  matrix  for  scattering  (S),  S,  can  be  setup 
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(15) 


where 


a,  =  Stokes  Vector  for  scattered  light 
cTj  =  Stokes  Vector  for  incident  light 


If  the  incident  light  is  right-circularly  polarized,  then  the  irradiancc  of  the  the  scattered  light  is 
related  to  Sn  +  Sj^  Also,  the  irradiance  Ij,  of  the  scattered  light  is  related  to  S^  -  S,^  The  Mueller 
scattering  matrix  element  S,^  can  be  expressed  as 

5i4=(/n-/L)/2/.  (16) 

where  I,  is  the  intensity  of  the  incident  light.  The  matrix  clement  S,^  is  readily  interpretable  in  terms 
of  the  difference  of  the  irradiances  of  scattered  light  for  incident  right-circularly  and  Icft-circularly 
polarized  light.  In  this  case  the  VCD  intensity  is  directly  proportional  to  the  S,^  matrix  element.  As  a 
result  a  VCD  theoretical  spectrum  can  be  used  to  predict  the  magnitude  of  the  Mueller  matrix  element 


3.  RESULTS 


The  Chemometric  and  Biometric  Modeling  Branch  (CBM)  of  CRDEC  has  undertaken  the  task  of 
generating  a  library  of  theoretically  determined  VCD  spectra  for  sugars,  sugar-like  molecules  and 
other  molecules  of  interest.  The  calculations  have  proceeded  in  the  following  stepwise  manner: 

1.  the  software  package  MMADS  (6)  was  run  using  the  option  MODEL  in  order  to  obtain  a 
crude  geometry 

2.  the  structure  obtained  in  step  (1)  was  then  optimized  using  MOPAC  version  3.0  (7) 

3.  optimized  structures  using  the  Gaussian  programs  Gaussian  82  (8-11)  or  CADPAC 
version  4.0  (12)  were  obtained.  (Both  packages  yield  identical  results  both  in  terms  of 
total  energy  and  optimized  structure.) 

4.  frequencies  and  rotational  strengths  were  obtained  by  running  CADPAC  version  4.0  (12) 
on  a  CRAY-XMP  computer. 

The  calculations  at  this  point  have  been  largely  carried  out  at  the  3-2 IG  level.  This  level  was  selected 
because  past  experience  indicates  very  reasonable  optimized  structures  are  found  as  compared  to 
experimental  structures.  As  is  well  known  with  ab  initio  determination  of  vibrational  motion,  the 
vibrational  frequencies  arc  typrcally  too  high  -  c.g.  the  C-H  stretch  predicted  for  3-21  G  level 
calculations  multiplied  by  a  scaling  factor  of  0.89,  found  empirically,  is  in  very  good  agreement  with 
experimental  data  (13). 

Attention  has  been  focussed  on  the  following  areas: 
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1.  determination  of  the  VCD  spectrum  for  molecules  for  which  literature  results  are 
available  -  the  molecules  selected  here  are  ethylene  oxide  (14),  the  isotopically 
substituted  ammonia  molecule  (NHDT)  (15,16),  and  methyl  thiirane  (17-19). 

2.  prediction  of  VCD  spectra  for  other  molecules  of  interest  -  hydrogen  peroxide,  HOOH, 
and  hydrogen  persulfide,  HSSH. 

3.  prediction  of  VCD  spectral  parameters  for  molecules  that  have  functionalities  similar  to 
sugars  (hydroxyacetaldehyde,  hydroxymethyl  ether,  ethylene  glycol  and  dimethyl  ether) 
and  sugars  with  3  carbon  atoms  (D-glyceraldehyde)  and  4  carbon  atoms  (D-erythrose, 
and  D-threose). 

Our  objectives  are  the  following: 

1.  develop  the  VCD  parameters  (frequency  and  rotational  strengths)  in  the  range  range 

800-1200  cm  *  that  indicate  transitions  of  high  rotational  strength  and  predict  the  correct 
sign. 

2.  develop  the  VCD  parameters  as  accurately  as  possible,  within  a  procedure  that  is 
computationally  reasonable,  for  sugars  in  order  that  the  results  can  be  dovetailed  into  the 
experimental  program  on  Mueller  Matrix  Spectroscopy  directed  by  one  of  the  authors 
(AHC). 


S.l  Ethylene  Oxide,  Isotopomer  of  Ammonia  (NHDT)  and  Methyl  Thiirane 

In  order  to  assess  the  feasibility  of  theoretically  predicting  VCD  spectra,  as  a  reasonable  first  step 
it  seemed  appropriate  to  carry  out  calculations  on  molecules  for  which  experimental  and  theoretical 
results  already  existed  in  the  literature.  First,  ethylene  oxide  has  been  previously  investigated  by 
Amos  et.  al.  (14)  The  geometry  used  by  these  investigators,  given  in  Figure  4,  was  used  here  also  and 
in  one  study  4-3 IG  level  calculations  were  performed  using  the  CADPAC  package.  The  atomic  polar 
tensor  and  the  atomic  axial  tensor,  given  by  equations  (4)  and  (5),  agree  with  Table  3  of  reference 
(14).  The  calculated  frequencies  and  rotational  strengths  determined  in  a  related  study  at  the  3-21  G 
level  are  reported  in  Table  2. 

Next,  the  isotopomer  of  ammonia  (NHDT)  was  studied  with  the  following  geometric  parameters: 
nitrogen  -  hydrogen  distance  (r^^  =  1.011609  A  and  hydrogen  -  nitrogen  -  hydrogen  bond  angle 
(^hnh)  =  106.697815  a.u.  Calculations  were  carried  out  at  both  the  3-2 IG  and  4-31  G  levels.  The 
results  are  reported  in  Tables  3  and  4.  The  energies  and  dipole  moments  agree  with  quantities 
reported  in  reference  (15).  An  important  point  to  note  here  is  that,  for  the  most  part,  signs  of 
transitions  with  high  rotational  strengths  remain  unchanged  with  basis  set  and  even  though  the 
magnitude  of  the  rotational  strength  changes  for  transitions  of  high  magnitude,  the  magnitude  remains 
high  as  the  basis  set  is  changed.  For  transitions  of  low  magnitude  the  sign  of  the  rotational  strength 
can  change  as  the  basis  set  is  changed.  The  frequencies  calculated  here  have  not  been  scaled  or 
adjusted  in  any  way  at  this  time.  This  scaling  of  the  force  field  is  discussed  as  an  important  point  for 
further  work  in  section  IV. 

Then,  attention  was  turned  to  methyl  thiirane.  Experimental  and  theoretical  work  was  carried  out 
by  Polavarapu  et.  al.  (17)  and  calculations  were  carried  out  by  Dothe,  Lowe  and  Alper  (18,19).  A 
comparison  of  results  are  given  Table  5  with  the  3-2 IG  optimized  geometry  given  in  Table  6  and 
shown  in  Figure  5.  Calculated  frequencies  are  qualitatively  similar  to  the  6-3 IG  level  calculation.  A 
comparison  of  the  rotational  strengths  with  the  rotational  strengths  determined  by  a  scaled  quantum 
mechanical(SQM)  force  field  determined  from  a  6-3 IG*  level  calculation  shows  the  agreement  is 
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reasonable  as  to  the  relative  magnitude  and  the  sign.  For  the  most  part  a  transition  corresponding  to  a 
high  rotational  strength  at  the  3-2 IG  level  of  calculation  agrees  in  sign  with  the  more  extensive  scaled 
6-3 IG*  level  of  calculation. 


S.2  Hydrogen  Peroxide  and  Hydrogen  Persulfide 

As  a  next  step  in  study  a  few  small  molecules  were  investigated.  Hydrogen  peroxide  (HOOH)  and 
hydrogen  persulfide  (HSSH)  were  investigated  with  the  results  reported  in  Tables  7-10  and  Figures  6 
and  7.  At  this  time  we  are  not  aware  of  any  experimental  VCD  spectra  of  HOOH  or  HSSH. 

S.S  Sugar-like  Molecules  and  Sugars 

Next,  attention  was  focussed  on  molecules  that  are  sugar-like  in  terms  of  functional  groups  and 
sugars  with  three  to  six  carbon  atoms  (19).  The  sugar-like  molecules  studied  were  isotopomers  of 
hydroxyacetaldehyde,  hydroxymethyl  ether,  ethylene  glycol  and  dimethyl  ether.  Calculated 
frequencies,  rotational  strengths  and  optimized  geometries  are  reported  in  Tables  11-18  and  the 
optimized  geometries  are  shown  in  Figures  8-11.  In  addition  optimized  structures  for  the  three  carbon 
atom  sugar  D-glyceraldehyde  and  the  four  carbon  atom  sugar  D-erythrose  have  been  obtained  and  are 
reported  in  Tables  19  and  20  and  shown  in  Figures  12  and  13.  The  frequencies  and  the  rotational 
strengths  have  not  been  reported  because  of  the  unavailability  at  this  time  of  the  CRAY-XMP 
computer  on  which  the  C  A  DP  AC  computer  program  is  run.  The  VCD  parameters  will  be  determined 
once  CRAY-XMP  computer  becomes  available  so  as  to  complete  this  phase  of  the  project.  In 
addition  optimization  of  the  geometry  of  the  four  carbon  atom  sugar  D-threose  is  nearing  completion 
and  once  this  is  completed  the  VCD  spectral  parameters  will  be  determined  from  CADPAC. 

S.^  Discussion 

Through  examination  of  given  spectral  ranges,  it  is  possible  to  scan  each  table  of  frequencies  and 
rotational  strengths  to  predict  at  which  frequency  a  maximum  rotational  strength  should  be  observed. 
The  quality  of  the  results  can  be  seen  by  focusing  on  the  VCD  results  found  for  methyl  thiirane  which 
arc  summarized  in  Table  5.  The  key  points  to  note  are  the  following: 

1.  The  calculation  is  sensitive,  because  of  the  nature  of  the  quantum  mechanical  operators 
involved  -  namely  the  angular  momentum,  to  the  choice  of  origin  of  the  coordinate 
system.  Comparisons  in  the  literature  at  present  find  better  agreement  with  the  choice  of 
origin  as  the  center  of  mass  of  the  molecule.  Very  key  to  the  calculation  is  the  level  of 
the  calculation.  As  mentioned  our  calculations  are  carried  out  at  the  3-2 IG  level,  the 
origin  chosen  was  the  center  of  mass. 

2.  The  inaccuracies  that  result  will  be  manifest  in 

•  force  field  generated 

•  atomic  polar  tensor 

•  atomic  axial  tensor 

Our  calculation  of  rotational  strengths,  compared  to  6-3 IG  level  calculations,  for  the 
majority  of  fundamental  frequencies  predict  transitions  of  high  rotational  strengths  and 
the  correct  sign.  Transitions  with  low  strengths  can  have  their  magnitude  and  sign 
changed  as  one  uses  another  basis  set. 
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The  calculations  on  ammonia  agree  with  results  reported  in  the  literature  as  to  the  values  of  the 
atomic  polar  and  atomic  axial  tensors  and  the  rotational  strengths.  The  frequencies  determined  in  the 
references  (15,16)  have  been  arrived  at  by  an  empirical  force  field  (Duncan-Mills).  As  has  been 
mentioned  the  ab  initio  unsealed  frequencies  will  be  somewhat  high  compared  to  experimental 
frequencies. 


4.  FUTURE  DIRECTIONS 

The  following  work  is  suggested  to  continue  the  project  toward  the  generation  of  a  library  of 
theoretically  predicted  VCD  spectra: 

1.  Complete  work  on  the  four  carbon  atom  sugars 

2.  Label  the  modes  of  vibration  for  all  molecules  that  were  studied  and  compare  with 
experimental  data  or  other  calculations  where  available. 

3.  Carry  out  studies  of  sugars  with  five  carbon  atoms  and  six  carbon  atoms.  Investigate 
sugars  with  ring  structures,  such  as  glucose. 

4.  Investigate  the  sensitivity  of  the  3*2  IG  level  of  calculation  to  carrying  out  optimizations 
with  greater  tolerances  for  gradient  convergence,  such  as  10-*.  Tlie  tolerance  used  for 
the  gradient  convergence  for  the  studies  reported  here  was  10-*. 

3.  Investigate  larger  basis  set  calculations,  such  as  6-31G  or  6*31G‘,  using  the  3-21G 
optimized  structures  as  starting  geometries. 

6.  Investigate  the  use  of  scaled  force  field  and  their  effects  on  calculated  frequencies, 
atomic  polar  tensors  and  atomic  axial  tensors  and  their  subsequent  efiect  on  rotational 
strengths. 


In  order  to  obtain  quantitative  accuracy  in  the  prediction  of  die  VCD  spectrum  of  a  molecule, 
some  form  of  scaling  appears  to  be  needed.  However,  if  the  main  intent  of  the  calculations  is  to 
predict  the  transitions  of  high  rotational  strength  and  their  correct  signs  then  3-2 IG  seems  to  be  an 
adequate  level  to  do  the  calculations  at  based  on  the  small  set  of  molecules  that  have  been  studied  up 
to  this  point  in  time. 


5.  SUMMARY 

We  have  undertaken  the  generation  of  a  library  of  theoretically  predicted  VCD  spectra.  The 
approach  involves  the  use  of  Stephens’  theoretical  formulation  of  the  rotation  J  strength  and  the  use  of 
Gaussian  calculations  at  the  3-2 IG  level.  The  results  on  molecules  studied  up  to  this  point  indicate 
that  for  transitions  corresponding  to  high  rotational  strengths  the  signs  remain  insensitive  to  the  basis 
set  level  used;  in  addition  the  magnitude  of  the  rotational  strength  of  a  transition  is  sensitive  to  the 
basis  set  level  that  is  used  but  a  transition  of  high  rotational  strength  at  one  level  of  calculation  will 
tend  also  to  be  high  at  another  level  of  calculation.  It  should  be  emphasized  that  the  key  assumptions 
being  used  are: 


1.  use  of  Stephens’  theoretical  formulation 


2.  use  of  Gaussian  3-2 IG  level  of  calculation 

The  level  of  calculation  affects  the  frequencies,  atomic  polar  tensor,  atomic  axial  tensor,  and  thereby 
the  rotational  strength. 


I 

I 

I 


-16- 


0.8  H 


■17' 


(R)-inethylthiirane  in  the  1500-800  cm-‘  region.  The  measurements 
were  made  with  neat  liquid  in  a  variable  path  length  cell  using 
minimum  path  lengths  15  /i).  This  figure  is  taken  from  the  reference  P. 
L.  Polavarapu  ct.  al.,  J.  Chem.  Phys.  86,  1 140-46  (1987).' 


Figure  2.  Schematic  diagram  of  an  experiment  in  which  radiation  is  transmitted  through 
an  optical  element.  The  detector  measures  the  irradiance  of  the  beam 
transmitted  by  the  polarizer  P  [see  C.  F.  Bohren  and  D.  F.  Huffman, 
'Absorption  and  Scatttering  of  Light  by  Small  Particles',  Wiley- 
Interscience,  New  York  (1983),  Chapter  2].' 


Figure  3.  Schematic  diagram  of  an  experiment  in  which  radiation  is  scattered  by  a 
single  particle  or  collection  of  particles.  The  collimated  detector 
responds  only  to  scattered  light  [see  C.  F.  Bohren  and  D.  F.  Huffman, 
'Absorption  and  Scattering  of  Light  by  Small  Particles',  Wiley- 
Interscience,  New  York  (1983),  Chapter  3]. 
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Figure  4.  The  geometry  of  ethylene  oxide. 

The  internal  coordinates  are  r^,,  =  1.086  A,  r^c  =  1,462  A,  r^^  =1,428 
^ocH  =  *  1^-5  ^ocH  =  ^  *»  ^HCH  =  ^  The  Cartesian  axes  are 

as  shown;  z  is  the  C,  molecular  axis  and  the  ring  is  in  the  yz  plane. 
The  origin  o  is  on  the  C,  axis,  0,853  A  from  the  oxygen  atom.  O  is  at 
the  center  of  mass  [  see  R.  D.  Amos,  N.  C.  Handy,  K.  J.  Jalkanen  and 
P.  J.  Stephens,  Chem.  Phys.  Letters  133 , 21-26  (1987)]. 


H6 


Figure  5.  Optimized  gcometry(3-21G  level)  of  methyl  thiirane  (CjHJS) 
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Figure  6.  Optimized  geometry(3-21G  level)  of  hydrogen  peroxide 


Figure  7.  Optimized  geometry(3-21G  level)  of  hydrogen  persulfide  (H^,) 
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Figure  11.  Optimized  gcometry(3-2!G  level)  of  dimethyl  ether  (CjH,0) 
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TABLE  1.  Stokes  Parameters  for  Polarized  Light 
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TABLE  2.  Calculated  Wavenumbers  v  and  Rotational 
Strengths  R  for  Ethylene  Oxide  (C^HP)  Based  on  Geometry  shown  in 
Figure  1.* 


V  (cm-‘) 

R  (x  10-**esu%:m’) 

935.32 

0. 

939.52 

0.000834 

972.85 

-0.000744 

1190.51 

0. 

1233.74 

0.001365 

1310.86 

0. 

1313.93 

-0.000218 

1316.13 

-0.000482 

1406.52 

0. 

1674.80 

0. 

1706.69 

0. 

3156.53 

0. 

3171.58 

0. 

3248.94 

0. 

3266.77 

0. 
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TABLE  3.  Calculated  Wavenumbers  u  and  Rotational 
Strengths  R  for  the  Isotopomer  of  Ammonia  (NHDT)  Based  on 
Geometry  (r„„  -  1.011609  A,  -  106.69781S  *)  used  by  Jalkanen 
et.  al.,  J.  Chem.  Phys.  142,  153(1987). 

[  The  total  energy  -  -55.870479593  a.u.  and  dipole  moment  = 
-0.855236  a.u.  ].  llie  Frequencies  and  Rotational  Strengths  Obtained 
by  Jalkanen  et.  al..  Based  on  a  Duncan-Mills  Force  Field,  are  Included 
in  Parentheses.  Calculations  are  Done  at  the  3-21  G  Level. 


v  (cm-‘)  R  (x  10-**  esu^m*) 


976.57(  826.62) 

1352.70(1195.18) 

1676.81(1483.96) 

2224.70(2176.42) 

2662.35(2600.81) 

3647.66(3561.51) 


-10.620(-9.29) 

19.432(18.99) 

-10.510(-9.70) 

0.090(0.00) 

0.787(0.43) 

-0.540(-0.19) 


TABLE  4.  Calculated  Wavenumbers  v  and  Rotational 
Strengths  R  for  the  Isotopomer  of  Ammonia  (NHDT)  Based  on 
Geometry  (r^^  =  1.011609  A,  *  106.697815  *)  used  by  Jalkanen 

et.  al.,  J.  Chem.  Phys.  142,  153(1987). 

(  The  total  energy  »  -56.102481085  a.u.  and  dipole  moment  = 
-0.904429  a.u.  ].  The  frequencies  and  rotational  strengths  obtained  by 
Jalkanen  et.  al.,  based  on  a  Duncan-Mills  force  field,  are  included  in 
parentheses.  Calculations  are  done  at  the  4-31  G  level. 


V  (cm-‘)  R  (x  10  **  esu^cm*) 


1022.67(  826.62) 
1357.39(1195.18) 
1677.65(1483.96) 
2219.84(2176.42) 
2658.03(2600.81) 
3641.00(3561.51) 


-12.327(-9.70) 
21.868(19.86) 
-11. 501  (-10.24) 
0.890(0.58) 
-0.319(-0.18) 
-0.060(-0.06) 
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TABLE  5.  A  Comparison  of  Wavenumbers  t/and  Rotational  Strengths  R  for  Methyl  Thiirane 


V  (cm-*) 

R  (x  10- 

esu^m’) 

Ref.  a 

Ref.b 

This  work 

Ref.  a 

This  work 

220 

238.23 

7.9 

2.196 

301 

305 

312.65 

-0.82 

-8.078 

377 

391 

427.46 

-1.7 

3.568 

585 

592 

538.24 

-14.4 

-9.637 

628 

644 

602.72 

22.1 

11.396 

861 

885 

925.39 

-36.8 

-14.875 

896 

957 

943.69 

54.2 

21.076 

908 

1011 

1043.99 

-21.6 

-19.318 

992 

1098 

1 108.86 

73.2 

21.285 

1043 

1151 

1172.80 

40.9 

52.519 

1078 

1197 

1207.44 

-173.7 

-98.711 

1151 

1224 

1248.83 

5.2 

-1.470 

1169 

1268 

1324.05 

35.9 

-11.557 

1341 

1446 

1506.98 

15.0 

46.377 

1387 

1543 

1586.01 

0.85 

-2.889 

1434 

1578 

1624.52 

17.7 

4.697 

1446 

1605 

1655.00 

-5.4 

8.194 

1452 

1614 

1664.66 

-3.5 

-14.743 

2861 

3089 

3205.64 

0.59 

-1.512 

2919 

3146 

3264.62 

-2.0 

7.475 

2943 

3168 

3272.64 

9.4 

-6.989 

3018 

3235 

3344.02 

3.6 

-0.427 

3040 

3276 

3378.51 

-10.5 

6.182 

3098 

3317 

3444.78 

3.2 

0.016 

Ref.  a  -  H.  Dothe,  MA.Lowe  and  J.S.Alper,  J.  Phys.  Ghent  92,  6246*49  ( 1 988). 

These  authors  have  used  the  scaled  quantum  mechanical  (SQM)  force 
field  to  obtain  vibrational  frequencies  and  the  transformation  matrix 
from  cartesian  coordinates  to  normal  coordinates  that  is  needed  in  the 
VCD  calculation  of  rotational  strengths.  Original  calculations  were 
carried  out  at  the  6-31  G*  level. 

Ref.  b  ■  P.  L.  Polavarapu  et.  al.,  J.  Ghent  Phys.  86,  1 140-46  (1987).  These 

authors  carry  out  ab  initio  calculations  at  the  6-31  G  level. 

This  work  -  Calculations  carried  out  at  the  3-21  G  level. 
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TABLE  6.  Optimized  Geometry  for  Methyl  Thiirane 
(C^HjCHg)  Based  on  3-2 IG  Level  of  Calculation 
[  The  total  energy  »  -511.9828737  a.u.  and  the  dipole  moment  = 
3.01852  D.  ] 


geometrical  coordinate  value  of  geometrical  coordinate 

r(Sj-C^ 

1.9345  A 

r(C,-C^ 

1.4636  A 

r(C,-HJ 

1.0708  A 

r(C,-HJ 

1.0700  A 

r(C,-HJ 

1.0715  A 

r(C,-C,) 

1.5082  A 

r(C,-H^ 

1.0835  A 

^C,-C,-S,) 

67.754  • 

118.408  * 

119.065  * 

116.724  * 

^C,-C,-C^ 

119.880  • 

^H,-C,-C,) 

110.412  * 

7(H,-C,-C3-S.) 

104.653  * 

-103.887  • 

i(H,-C,-C,-S.) 

103.046  • 

-107.542  • 

t(H,-C,-C,-S^ 

76.805  • 

t(H,-C,-C,-Sj) 

-163.473  • 

i(H„-C,-C,-S^ 

-42.922  • 
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TABLE  7.  Calculated  Wavenumbers  v  and  Rotational 
Strengths  R  for  Hydrogen  Peroxide  (HOOH)  Based  on  3-2 IG  Level 
Optimized  Geometry 


mode 

0  (cm-‘) 

R  (x  10  “esu\:m*) 

torsion 

87.10 

-6.117 

1 167.99 

-0.174 

1278.05 

-0.506 

1644.43 

0.094 

3328.40 

-1.236 

3845.25 

1.888 

TABLES.  Optimized  Geometry  for  Hydrogen  Peroxide  (HOOH)  Based  on  3-2 IG  Level  of 
Calculation 


geometrical  coordinate 

value  of  geometrical  coordinate 

r(O-H) 

0.9703  A 

r(O-O) 

1.4731  A 

^H-O-O) 

99.451  • 

»(H-0-0-H) 

178.595  * 
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TABLE  9.  Calculated  Wavenumbers  v  and  Rotational 
Strengths  R  for  Hydrogen  Persulfide  (HSSH)  Based  on  3*210  Level 
Optimized  Geometry 


mode 

V  (cm-*) 

R  (x  10-“csu^cm*) 

torsion 

337.34 

-159.66 

462.28 

-0.808 

897.81 

30.828 

912.98 

■29.138 

2632.15 

-6.583 

2633.82 

0.224 

TABLE  10.  Optimized  Geometry  for  Hydrogen  Persulfide  (HSSH)  Based  on  3*2  IG  Level  of 
Calculation 


geometrical  coordinate 

value  of  geometrical  coordinate 

r(S-H) 

1.3514  A 

r(S-S) 

2.2638  A 

^H-S-S) 

96.707  * 

»(H-S-S-H) 

93.830  * 
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TABLE  11.  Calculated  Wavenumbers  t/and  Rotational  Strengths  R  for 
the  Isotopomer  of  Hydroxyacetaldehyde  (CHDOH>CHO)  Based  on  3- 
21 G  Level  Optimized  Geometry 


1/  (cm-‘)  R  (x  10-**  esu^m*) 

174.52 

-1.5369 

317.00 

-0.8014 

397.23 

-9.3339 

755.26 

67.952 

808.78 

-28.285 

882.65 

-81.512 

1004.03 

27.579 

1166.53 

74.788 

1222.99 

-57.161 

1389.71 

0.3105 

1487.43 

43.712 

1519.10 

-42.697 

1541.64 

-5.9977 

1906.71 

10.576 

2358.19 

24.699 

3206.11 

-63.704 

3223.24 

39.647 

3828.49 

1.2695 
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TABLE  12.  Optimized  Geometry  for  the  Isotopomer  of  Hydroxyacetaldehyde 
(CHDOH-CHO)  Based  on  3-2 IG  Level  of  Calculation. 

[  The  total  energy  =  -226.484421186  a.u.  and  the  dipole  moment  = 
2.9617  D.  ] 


geometrical  coordinate  value  of  geometrical  coordinate 

r(0,-C^ 

1.4255  A 

r(C,-C,) 

1.5094  A 

r(0,-H^ 

0.9701  A 

r(C,-HJ 

1.0849  A 

r(C,-HJ 

1.0849  A 

r(C,-0,) 

1.2103  A 

r(C,-H^ 

1.0826  A 

<KC,-C,-0,) 

110.420  • 

108.551  • 

^H,-C,-0.) 

110.663  • 

110.669.' 

^0,-C.-CJ 

121.590  * 

116.115  * 

»(H,-C,-0.-C,) 

-0.025  * 

7(H,-C,-0,-H^ 

120.293  • 

-120.360  • 

i(0,-C,-C,-0^ 

0.003  • 

i(H,-C,-C,-0^ 

179.997  * 
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TABLE  13.  Calculated  Wavenumbers  u  and  Rotational  Strengths  R  for 
the  Isotopomer  of  Hydroxymethyl  ether  (CH,-0-CHDOH)  Based  on 
3-2 IG  Level  Optimized  Geometry 


V  (cm-‘)  R  (x  10-**esu^m*) 

150.63 

25.440 

197.25 

-3.754 

136.10 

17.703 

419.20 

-33.412 

618.97 

17.865 

924.28 

19.796 

991.27 

74.355 

1072.03 

22.944 

1150.26 

-106.04 

1224.07 

38.492 

1268.74 

18.074 

1313.02 

20.976 

1411.83 

-50.363 

1524.23 

6.133 

1551.44 

-63.669 

1625.46 

8.489 

1671.00 

-6.482 

1698.14 

8.946 

2398.46 

-10.482 

3209.41 

8.844 

3263.13 

5.903 

3271.06 

-31.401 

3285.84 

4.546 

3867.18 

24.137 
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TABLE  14.  Optimized  Geometry  for  the  Isotopomer  of  Hydroxymethyl  ether 
(CH,-0-CHD0H)  Based  on  3-2 IG  Level  of  Calculation. 

[  The  total  energy  »  -227.661 102  a.u.  and  the  dipole  moment  =  0.37208 
D.] 


geometrical  coordinate 

value  of  geometrical  coordinate 

r(C,-0^ 

1.4418  A 

r(0,-C,) 

1.4143  A 

r(0,-CJ 

1.4210  A 

r(Oi-CJ 

1.0814  A 

T(0,-CJ 

0.9671  A 

^c,-o,-c^ 

114.688  • 

^0,-C,-0^ 

112.470  • 

106.572  • 

110.816  * 

^H,-C.-0,) 

110.830  • 

111.513  * 

105.479  * 

<K»,oO,-C^ 

110.488  • 

i(0,-C,-0,-C,) 

62.680  • 

-178.097  • 

i(H.-C,-0,-C,) 

62.801  • 

i(H,-C,-0,-C,) 

-58.252  • 

f(H,-C,-0,-C,) 

-55.079  • 

i(H,-C,-0,-C^ 

-174.840  * 

7(H,„-0,-C,-0^ 

60.925  • 
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TABLE  15.  Calculated  Wavenumbers  v  and  Rotational  Strengths  R  for 
the  Isotopomer  of  Ethylene  Glycol  (CHDOH-CH,OH)  Based  on  3-2 IG 
Level  Optimized  Geometry 


j/ (cm-‘)  R  (x  10-**esu^m*) 

145.07 

13.598 

275.29 

220.23 

287.60 

140.72 

294.46 

-384.63 

488.68 

6.512 

810.61 

-9.441 

979.33 

24.161 

1043.30 

26.607 

1128.76 

4.467 

1144.91 

-36.410 

1232.74 

29.689 

1323.41 

-12.872 

1335.22 

-3.809 

1397.95 

-28.120 

1489.16 

31.741 

1533.26 

-24.860 

1581.46 

-3.362 

1704.85 

1.119 

2396.71 

-6.257 

3212.47 

4.063 

3234.02 

51.990 

3258.77 

-55.673 

3872.00 

6.592 

3877.91 

-3.438 
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TABLE  16.  Optimized  Geometry  for  the  Isotopomer  of  Ethylene  Glycol 
(CHDOH-GHjOH)  Based  on  3-21G  Level  of  Calculation. 

[  The  total  energy  -  -227.6518901  a.u.  and  the  dipole  moment 
2.56172  D.  1 


geometrical  coordinate 

value  of  geometrical  coordinate 

r(0,-C^ 

1.440  A 

r(C,-C,) 

1.521  A 

r(C,-OJ 

1.440  A 

r(0,-H^ 

0.966  A 

r(C,-H^ 

1.082  A 

^C,-C,-0^ 

110.279  • 

<i(0,-C,-C,) 

105.448  • 

110.994  * 

109.407  • 

^H,-C,-C,) 

109.407  • 

109.407  * 

109.407  • 

110.994  * 

i(0,-C,-C,-OJ 

178.917  • 

i(H,-0,-C,-C,) 

78.201  •• 

i(H,-C,-C,-0^ 

-56.813  • 

KH,-C,-C.-0^ 

61.786  • 

57.886  • 

»(H,-C,-C,-0^ 

-60.910  • 

-176.349  • 
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TABLE  17.  Calculated  Wavenumbers  »/and  Rotational  Strengths  R  for 
the  Isotopomer  of  Dimethyl  Ether  (CHOT-O-CH,)  Based  on  3-2 IG 
Level  Optimized  Geometry 


v(cm-‘)  R  (x  10  **csu^m*) 

165.27 

-5.2052 

258.97 

-7.8217 

393.57 

7.6908 

883.63 

-3.7270 

956.41 

-11.465 

1055.51 

9.7607 

1141.22 

0.4384 

1236.85 

-5.4182 

1280.68 

-0.2552 

1305.47 

-2.9711 

1447.48 

26.488 

1486.33 

-19.965 

1635.26 

1.1278 

1678.55 

1.0211 

1690.78 

0.6431 

1989.75 

13.079 

2360.74 

-14.611 

3181.62 

0.2114 

3231.41 

0.1558 

3238.94 

-1.1410 

3252.97 

0.2455 
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TABLE  18.  Optimized  Geometry  for  the  Isotopomer  of  Dimethyl  Ether 
(CHDT-O-CHg)  Based  on  3-2 IG  Level  of  Calculation. 

[  The  total  energy  -  -153.212180696  a.u.ole  moment  ■ 

1.8925  D.  ] 


geometrical  coordinate  value  of  geometrical  coordinate 

r(C,-0^ 

1.4324  A 

r(0,-CJ 

1.4324  A 

r(C,-HJ 

1.0840  A 

r(C.-H^ 

1.0840  A 

r(C,-H^ 

1.0840  A 

r(C,-H,) 

1.0840  A 

r(C,-H^ 

1.0840  A 

r(C,-H^ 

1.0840  A 

^C,-0,-Cj) 

113.190* 

110.001  * 

110.001  * 

110.001  * 

110.005  * 

1 10.005  * 

110.005  * 

»(H,-C.-0,-C,) 

-179.981  * 

i(H,-C,-0,-C,) 

59.921  * 

t(H.-C.-0,-C,) 

-59.884  * 

i(H,-C,-0,-0^ 

179.993  * 

i(H.-C,-0,-C,) 

59.897  * 

i(H,-C,-0,-CJ 

-59.910  * 
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TABLE  19.  Optimized  Geometry  for  the  Isotopomer  of  D-Glyceraldehyde 
(CHO-HCOH-GHjOH)  Based  on  3-2 IG  Level  of  Calculation. 

[  The  total  energy  =  -339.7416503  a.u.  and  the  dipole  moment  =  3.4498 


geometrical  coordinate  value  of  geometrical  coordinate 

r(C,-CJ 

1.5103  A 

r(C,-C^ 

1.5286  A 

r(C.-HJ 

1.0820  A 

r(C,-OJ 

1,2077  A 

r(C,-HJ 

1.0810  A 

r(C,-0,) 

1.4439  A 

r(C,-H0 

0.9655  A 

r(C3-Hj 

1.0827  A 

r(C3-H,) 

1.0779  A 

r(C,-0,,) 

1.43441  X 

r(Ou-HJ 

0.9681  A 

^C,-C,-C») 

109.217  • 

113.171  • 

<KO,-Ccc;) 

123.679  * 

1 10.335  • 

^o,-c,-c,) 

109.089  • 

^H,-0,-C,) 

112.428  • 

^H,-C.-C,) 

108.690  • 

110.831  • 

<i(0,,-Cj-C,) 

108.504  • 

^H„-0„-C,) 

107.757  • 

i(H,-C,-C,-C,) 

78.170  • 

-99.786  • 

i(H,-C,-C^-HJ 

200.808  • 

-33.907  • 

7(H,-0,-C,.C,) 

-73.590  • 

<H,-C,-C,-C.) 

-181.367  • 

i(H.o-C,-C,-C.) 

58.027  • 

HO,,-C,-C,-C,) 

-59.027  • 

-59.687  ' 

TABLE  20.  Optimized  Geometry  for  D-Erythrose  (CHO-HCOH-HCOH-CH,OH)  Based  on 
3-2 IG  Level  of  Calculation. 

[  The  total  energy  =  -453.008065  a.u.  and  the  dipole  moment  «  3.0047 1 
D.] 


geometrical  coordinate 

value  of  geometrical 

r(C.-C^ 

1.5156  A 

r(C,-C,) 

1.5215  A 

r(C,-C; 

1.5226  A 

r(C^-HJ 

1.0839  A 

TiC,-0^ 

1.2074  A 

r(C,.H,) 

1.0827  A 

r(C,-0; 

1.4282  A 

r(0,-HJ 

0.9714  A 

r(C,-HJ 

1.0807  A 

r(C,-0„) 

1.4438  A 

r(Ou*H.,) 

0.9698  A 

r(C,-HJ 

1.0821  A 

r(C,-HJ 

1.0771  A 

r(C,-OJ 

1.4509  A 

0.9663  A 

^C,-C,-CJ 

110.196  • 

^C,-C,-C^ 

112.421  • 

113.564* 

^0,-C,-C^ 

123.165  * 

110.846* 

^0,-C,-C,) 

109.653  * 

106.172  * 

109.921  * 

^o„-c,-c,) 

104.007  * 

107.836  * 

111.752  * 

^o,,-c,-c,) 

108.925  * 

^o».-c,-c.) 

107.857  * 

111.275  * 

t(C,-C,-C,-C,) 

171.623  * 

t(H.-C,-C,-C,) 

137.045  * 

1(0,-C,-C,-C^ 

-43.836  * 

i(H,-C,-C,-CJ 

51.303  * 

i(0,-C,-C,-CJ 

-67.691  * 

81.697  * 

46.310  * 

i(0„-C,-C,-C  j 

-72.685  * 

i(H„-0.,-C.-C,) 

203.203  * 

-41- 


HH„-C«-C.-C^ 

-73.658  * 

7(H„-C,-C,-C^ 

47.640  * 

1(0.,-C,-C,-C,) 

-197.374  • 

t(H  ,,-0  ij"C^-C 

86.553  • 

-42- 
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